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{NJ ■ Abstract 



In this talk it is shown how nonrelativistic QCD (NRQCD) can be used to determine next-to-next- 



to-leading order relativistic and short-distance contributions to the total tt production cross section 
in the threshold regime at lepton colliders. A recipe for the calculation of all such contributions 



for the total photon mediated production cross section is presented and a review of the already 
known Abelian next-to-next-to-leading results is given. 



X. 

■ Introduction 

The production of tt pairs in the threshold region at future lepton colliders like the NLC (Next Linear 
Collider) or the FMC (First Muon Collider) offers a unique opportunity to carry out precision tests 
of QCD in a completely new environment. Due to the large top mas^] [Mt ~ 175 GeV), which allows 
for the decay channel t — > Wb, hadronization effects can be neglected in a first approximation 
This makes the ti production cross section in the threshold regime (including various distributions) 
calculable from perturbative QCD (and electroweak interactions), which then allows for precise ex- 
tractions of the top quark mass and the strong coupling once the cross section is measured. In fact, 
experimental simulations for the NLC and the FMC |2|, [3| have shown that experimental errors of 
around 100 MeV for the top quark mass and of around 0.002 for a s (M z ) can be expected for a cross 
section measurement with a total integrated luminosity of 50 — 100/6 -1 . In particular the prospect 

'Invited talk presented at the Workshop on Physics at the First Muon Collider and the Front End of a Muon Collider, 
Fermilab, November 6 - 9, 1997, to appear in the proceedings. 

1 Throughout this talk M t is understood as the top quark pole mass. 



for the error in the top quark mass measurement beats any hadron collider experiment. However, 
the errors for Mt and a s given above do not contain any theoretical uncertainties. At this point it 
is illustrative to recall that the standard present day formalism used for describing ti production at 
threshold consists of solving the nonrelativistic Schrodinger equation with a QCD potential which for 
small distances is given by perturbative QCD up to one loop ||], |j| and for large and intermediate 
distances by fits to quarkonia spectra (and leptonic decay widths) || [7j . The results are then modified 
by various 0(a s ) short-distance corrections which makes the results correct at the next-to- leading 
order (NLO) level, i.e., they properly include all 0{a s ) corrections.^ NNLO {0{a 2 s )) corrections have 
never been taken into account so far. Their contributions, however, can be sizable. As an example, 
consider the 0{a 2 s ) relativistic corrections to the total cross sections which can lead to a shift in the 
location of the IS peak of order M t a^ ~ 150 MeV and a corrections of order a 2 ~ 3% in the size 
of the cross section (for a s ~ a s (M t a s ) ~ 0.17). Further, even the 0(a 2 ) short-distance corrections 
normalizing the total cross sections might be large if the huge size of the 0{a s ) corrections of order 
—20% is taken into account. From this point of view it is clear that the theoretical uncertainties in 
the present day analyses are certainly not negligible and that full control over all NNLO effects should 
be gained. 

Unfortunately the formalism described above is constructed in a way that makes a systematic 
and rigorous implementation of all NNLO effects from first principles QCD conceptually difficult if 
not impossible - a consequence of the use of phenomenological information in the potential for large 
and intermediate distances. In principle, this formalism has to be considered as a (sophisticated) 
potential model approach which cannot be improved in a rigorous way at all. I therefore propose to 
rely on perturbative QCD only. This means that one applies the perturbative QCD potential for all 
distances. In fact, such a decision seems to be just natural if one takes into account that the ti system 
is almost insensitive to the large distance (i.e. non-perturbative) contributions in the QCD potential. 
This makes the framework in which effects beyond the NLO level shall be determined more transparent 
and still leaves the possibility to incorporate the non-perturbative effects later as a perturbation (see 
e.g. Q for an approach of this sort). 

In this talk I demonstrate how NRQCD || can be used to determine NNLO relativistic and 
short-distance corrections to the total ti production cross section. For simplicity only the production 
through a virtual photon is considered. The generalization for production through different currents 
is straightforward. I want to stress that I do not talk about the peculiar NNLO finite width effects 
coming from the off-shellness of the top quark, time dilatation effects and the interaction among the 



decay products. These effects have been addressed previously in a number of publications JlOj, [l_l|, |12fl , 
but no rigorous and consistent description of them for the ti cross section has been found yet. To 
achieve full NNLO accuracy these finite width effects should eventually be included. For now they 
remain as an unsolved problem. As far as the NNLO relativistic corrections discussed in this talk are 
concerned I will use the naive replacement 

E = ^-2M t — >E = E + iT t (1) 

in the spirit of JlJ in order to examine their size and properties, where Tt represents a constant which 
is not necessarily the decay width of a free top quark. 

2 The solutions of the nonrelativistic Schrodinger equation with the one-loop corrected QCD potential contains, in 
the language of Feynman diagrams, the resummation of terms oc (a a /v) n x [l,a a ], n = 0, 1, . . . ,oo, (v being the top 
quark cm. velocity) to all orders in a s . Because in the threshold region \v\ < a 3 we count all terms oc a 3 /v of order one. 
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The tt Cross Section in NRQCD 



NRQCD is an effective field theory of QCD designed to handle nonrelativistic systems of heavy quark- 
antiquark pairs to in principle arbitrary precision. It is based on the separation of long- and short- 
distance effects^ by reformulating QCD in terms of a non-renormalizable Lagrangian containing all 
possible operators in accordance to the symmetries in the nonrelativistic limit. The NRQCD La- 
grangian reads 



-NRQCD 



q=u,d,s,c,b 



D 2 t a 3 g 
iDf + ai h ao q H cr ■ is 

2M t 8M ( 3 2M t 



8M? 1 ; 8M? K ' 



!& + .... (2) 



The gluonic and light quark degrees of freedom are described by the conventional relativistic La- 
grangian, whereas the top and antitop quarks are described by the Pauli spinors ip and Xi respectively. 
For convenience all color indices are suppressed. The straightforward antitop bilinears are omitted 
and only those terms relevant for the NNLO cross section are displayed. Dt and D are the time and 
space components of the gauge covariant derivative and E l = G 0t and B i = ^e l ^ k G^ k the electric 
and magnetic components of the gluon field strength tensor (in Coulomb gauge). The short-distance 
coefficients 01, . . . , 05 are normalized to one at the Born level. 

To formulate the normalized total ti production cross section (via a virtual photon) R = 
a(^ + ^_ — > 7* — > tt)/a p t (a p t = 47ro; 2 /3s) in the nonrelativistic region at NNLO in NRQCD we start 
from the fully covariant expression for the cross section 

R tf) = Im[ ( 1 T J M (g) (-g) ] >] , (3) 

where Qt = 2/3 is the electric charge of the top quark. We then expand the electromagnetic current 
(in momentum space) J„(rbg) = (ij^fy^q) which produces/annihilates a pair with cm. energy \fq 1 ~ 
in terms of 3 Si NRQCD currents up to dimension eight (i = 1, 2, 3)0 

ji(q) = 61 (^W)(«) " (^H Dfx)(q) + -.., (4) 

where the constants 61 and 62 are short-distance coefficients normalized to one at the Born level. 
Inserting expansion @ back into Eq. (@) leads to the NRQCD expression of the nonrelativistic cross 
section at the NNLO level 
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^nnlo(^) = ~T7T Cll/^hard, Mfac) I m <Al(E, // so ft, A*fa 

C2(Mhard,A*fac)Im Az{E, /isoft , A'fac) +•••, (5) 



3M t 4 



3 In this context "long-distance" is not equivalent to "non-perturbative" . 

4 Only the spatial components of the current contribute. The expansion of jtj,(—q) is obtained from Eq. (Q) via charge 
conjugation symmetry. 
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where 



Ay = (0|(^x) (x^^lO), (6) 

A 2 = ±(0\(ftffx) (x^(-iDN) +h.c.|0). (7) 

The cross section is expanded in terms of a sum of absorptive parts of nonrelativistic current-current 
correlators (containing long-distance physics) multiplied by short-distance coefficients C% {% = 1, 2, . . .). 
In Eq. (|D I have also shown the dependences on the various renormalization scales: the soft scale 
// ao f t and the hard scale /Uhard are governing the perturbative expansions of the correlators and the 
short-distance coefficients^], respectively, whereas the factorization scale /Uf ac essentially represents the 
boundary between hard (i.e. of order Mt) and soft momenta. Because this boundary can in principle 
be chosen freely, both, correlators and the short-distance coefficients, in general depend on it (leading 
to new anomalous dimensions). Because the term in the second line in Eq. (||) is already of NNLO (i.e. 
suppressed by v 2 ) we can set C2 = 1 and ignore the factorization scale dependence of the correlator 
Ai- The calculation of all terms in expression (||) proceeds in two basic steps. 

1. Calculation of the nonrelativistic correlators. - Determination of the correlators in Eq. (||) by 
taking into account the interactions up to NNLO displayed in £ N rqcd- 

2. Matching calculation. - Calculation of the constant C\ up to O(o%) by matching expression (||) 
to the fully covariant cross section at the two-loop level in the (formal) limit a s -C v <C 1 where 
an expansion in (first) a s and (then) v is feasible. 

Calculation of the nonrelativistic correlators: In Coulomb gauge the gluon propagation is separated 
into a longitudinal, instantaneous (i.e. energy-independent) and a transverse, non-instantaneous (i.e. 
energy-dependent) propagation. The longitudinal exchange between the tt pair is described by an 
instantaneous potential. (The Coulomb potential is just the LO interaction caused by the longi- 
tudinal exchange.) The transverse exchange, however, leads to a temporally retarded interaction, 
closely related to Lamb-shift type effects known in QED. Fortunately, because the tt pair is pro- 
duced/annihilated in a color singlet (S-wave) configuration, the energy dependence of the transverse 
gluon exchange leads only to NNNLO (i.e. O(a^) relative to the effects of the LO Coulomb exchange) 
contributions.^ We therefore ignore the energy dependence of the transverse gluons which allows us to 
formulate all interactions contained in the NRQCD Lagrangian in terms of instantaneous potentials. 
In other words, as far as the nonrelativistic correlators at NNLO in Eq. (^) are concerned, NRQCD 
reduces to a two-body (top-antitop) Schrodinger theory. The potential in the resulting Schrddinger 
equation is determined by considering tt — * ti one gluon exchange t-channel scattering amplitudes 

5 The scales /i so f t and /ihard arise from the light degrees of freedom in £nrqcd and are already present in full QCD, 
whereas ^if ac is generated by "new" UV divergences in NRQCD diagrams. It is crucial to consider fi so ft and /Xhard as 
independent scales. Because both, nonrelativistic correlators and short-distance coefficients, are calculated perturbatively 
a residual dependence on fi ao ft and /Xhard remains. 

6 This can be seen by either using formal counting rules (see e.g. ijlj] , (mJ) or from positronium results where this 
phenomenon is well known. From the physical point of view the suppression comes from the fact that the transverse 
gluon is radiated after the ti pair is produced and absorbed before the ti pair is annihilated. This process is already 
suppressed by v 2 due to the dipole matrix element for transverse gluon radiation/absorption. If the gluon also carries 
energy, another (phase space) factor v arises because the gluon essentially becomes real. For the same reason no top 
quark self energy or crossed ladder diagrams have to be considered. 
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in NRQCD. To NNLO (i.e. including potentials suppressed by at most a 2 , a s /M t or 1/M 2 relative 
to the Coulomb potential) all potentials are already known and read (a = Cf a s (/i so ft), Cf = 4/3, 
C A = N c = 3, r = \r\) 



air 



5^{r) + 



V 2 + 



1 



(fV)V 



3a r 1 



M 4 2 r 3 



+ 



3 a 



2 M 2 r 3 



L(S t + S { 



C F 2 Mi r 2 



(8) 
(9) 



where St and £f are the top and antitop spin operators and L is the angular momentum operator. V BF 
is the Breit-Fermi potential known from positronium and V NA a purely non-Abelian potential generated 
through non-analytic terms in one-loop NRQCD (or QCD) diagrams containing the triple gluon vertex 
(see e.g. |l5| for an older reference). The Coulomb potential V c (r) = —a/r [1+ corrections up to 0(a 2 )] 
is not displayed due to lack of space. Its 0(a s ) and 0(a 2 ) corrections have been determined in [||, |j 
and |[q] , respectively. 

The nonrelativistic correlators are directly related to the Green function of the Schrodinger 
equation 



V 



V 



TT ~ 777s + V c(7) + ^ F (r) + V NA (7) - E G(r, 7, E) 



M t AM? 



5®{ 



r — r 1 ) 



where V BF is evaluated for the 3 S\ configuration only. The correlators read 

Ax = 6N C \ lim G{f,r,E)}, 

L |r|,|r*|— »0 J 

A 2 = M t EAx. 



(10) 

(11) 
(12) 



Relation ( JTT] ) can be easily inferred by taking into account that the Green function G(r, r* ', E) describes 
the propagation of a top- antitop pair which is produced and annihilated at distances \f\ and |r"|, 
respectively. (A more formal derivation can be found e.g. in ||.) Relation (|i~2]), on the other hand, is 
obtained through the equation of motion (10). Because the exact solution of Eq. (|l0[ ) seems to be an 



impossible task and we start with the well-known Coulomb Green function G^ |17]] iVc' 1 {r) = —a/r) 



r(0). 



(13) 



and incorporate the higher order terms using Ray leigh- Schrodinger time-independent perturbation 
theory (TIPT). It should be noted that the limit |r|, \r'\ — > causes UV divergences which have to be 
regularized using either an explicit short-distance cutoff or dimensional regularization. The freedom in 
the choice of the regularization parameter is the origin of the factorization scale dependence mentioned 
before. 

Matching calculation: The determination of the short-distance coefficient C\ up to C(a 2 ) can be carried 
out in two ways: One either calculates the constants b\ in Eq. (Q) through two-loop matching at the 
amplitude level for the electromagnetic vertex in full QCD and NRQCD or one matches expression (||) 



5 



directly to the two- loop cross section calculated in full QCD. Both ways of matching must be carried 
out for stable top quarks (Tt = 0) and are performed in the (formal) limit a s <C v <C 1 where NRQCD 
and full QCD are applicable in the conventional multiloop approximation .[] In our case one has to 
match at the two-loop level including terms up to NNLO in an expansion in v. Technically the second 
way of matching, called "direct matching" IS], is simpler because it allows for a sloppier treatment 
of the regularization procedure. (In fact, using the first way one has to be very careful to use exactly 
the same regularization for the matching calculation as in the calculation of the correlators. This is 
a quite tricky task if one wants to avoid solving the Schrddinger equation in D dimensions.) The 
disadvantage of direct matching, however, is that matching is carried out at the level of the final 
result which practically eliminates the possibility to use the calculated short-distance coefficient for 
any other process. Further, it requires that a multiloop expression for the cross section is at hand. 
For convenience, I use the "direct matching" method. 



Some Explicit Results 



This talk would not be complete without making the somewhat general discussion before more explicit. 
In the following I will carry out the program described above for all Abelian contributions, i.e., for 
those effects which also exist in QED. I will not present any technical details and refer the interested 
reader to [p2Ql , where the calculations have actually been carried outj^j I start with the well-known 
expression for the Coulomb Green function IT/ 

r2 



2vr 



dte 



2iM t vrt 



1 + t 
t 



E 

Mt 



(14) 



The NNLO corrections coming from the kinetic term — V 4 /4M 3 and the Breit-Fermi potential V B 
are calculated via TIPT and read 



5G(0,0,E) 



r V 4 

dx 3 G { V(0,x,E) 



4M t 3 



V BF (x)] Gf\x,0,E). 



(15) 



Abelian corrections coming from the one- and two-loop contributions in the Coulomb potential do not 
exist because we define a s = o^j§ ■ As mentioned before, taking the limit \r\, | ■r 1 ' | — > in Eq. (14) 
and the integral (^) leads to UV divergences which I regularize with the short-distance cutoff //f ac - 
This leads to the following result for the correlator A\ at NNLO 



A! 



2vr 



i v 



ln(- 



-i— — )+7 + * 1 



+- 



9 M, 



2vr 



i [v H — v" 



/"fa 



a[ 1 + 2v 2 



2v 



(16) 



ln(- 



Mtv, 
-a— H + 7 + * 1 

Wac 



2v 



where 7 is the Euler constant and the digamma function. All power divergences oc /if ac /Mt are 
freely dropped and /if ac is defined in a way that expression ([jl]) takes the simple form shown above. 
For the correlator A2 only the LO contribution in (|i~6|) is relevant and we arrive at 



Ao 



9M t 4 
2vr 



1 v 



ln(- 



-t— M + 7 + * 1 



(17) 



Mfac ' V 

7 Other matching points like v = are also possible, see e.g. [ft8||). 

8 In Ref. |2(]] the calculations have not been formulated in the framework of NRQCD. The results, of course, are not 
affected by this. 
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There are no non-Abelian contributions to A%- The Abelian contributions to C\ are calculated by 
expanding expression (|5|), expanding (first) for small a s and (then) v, keeping terms up to order a 2 , 
and NNLO in the v expansion,^ and demanding equality to the corresponding two-loop cross section 
calculated in full QCD |(| || for T t = (v = (E/M t ) 1/2 ), 



pNNLO 
- fl 21oop QCD 



N c Ql 



ci 



3 

— v 
2 



7T 



17 
16 



+ 



TT 



3tt 2 



7T n 
6u + — V 



3Ci + 



49C£vr 
192" 



2 3 



Ci Lav 



where C A 
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1)]. The result for C\ then reads 



Ci 



1 



|r , «s(^hard) , o, s 
40f h a s (^hard) 



7T 



CAbcl + ~ C p 



In 



(-) 



(18) 



(19) 



Expression ( |i9[ ) does not contain any energy dependent (or even IR divergent) contributions because 
NRQCD and QCD have the same low energy behavior, i.e., all the energy dependence is contained in 
the correlators. Apart from the dependence on the hard scale /ihardj which remains because [Ci] Abcl 
represents a truncated perturbative series,^ there is a dependence on the factorization scale //f ac . As 
can be seen in Eq. (|i~6|), for T t = this dependence is cancelled by a corresponding /Zf ac -dependent 
term in [^4i]Abei-Q For Tt ^ there is a small contribution cx a^hx^- which is not cancelled. 
This ambiguity arises from our ignorance of a consistent treatment of the NNLO finite width effect 
mentioned at the beginning. In Fig. |l| the relative size (in %) of the NNLO Abelian contributions 
C^nnlo - -^nloV-^nlo (^nlo contains only the contribution from the Coulomb Green function, Eq. fll4|), 
in A\ and the terms in C\ up to 0(a s )) is plotted in the energy range —10 GeV< E < 10 GeV for 
M t = 175 GeV, a s {M z ) = 0.118 and r t = 1.56 GeV (solid line)/0.80 GeV (dashed line). For the scales 
the choices ^ so f t = 30 GeV and /^hard = Mfac = Mt have been made and two-loop running of the strong 
coupling has been used. It is evident that the corrections are indeed at the level of several percent 
and that they are fairly insensitive to the value of Tt indicating that the size of the Abelian NNLO 
contributions is not affected by the ignorance of a consistent treatment of the finite width effects. I 
also would like to note that the largest source of theoretical uncertainty in the cross section -R^nlo 
comes from the dependence on the soft scale /i so ft - clearly because it governs the perturbative series 
describing the long-distance (i.e. low energy) effects for which the convergence can be expected to be 
worse than for the short-distance coefficients. A thorough examination of this behavior, however, has 
to wait until all the NNLO relativistic corrections are calculated. 



At this point we set /Lt so ft = /ihard because for a, C s « 1 a distinction between the soft and the hard scale is 
irrelevant. 

10 The dependence on /ihard of the 0(a s ) term in Eq. ( p^| is not cancelled by terms in the 0(a%) contributions because 
non-Abelian and massless quark corrections are not considered. 

11 The invariance under changing the factorization scale ^t lac can be used to resum renormalization group logarithms. 
However, I would like to warn the reader to blindly apply renormalization group methods in the belief this would 
represent a resummation of "leading logarithms" . Although it is true that a naive resummation of logarithms is possible 
in this way, the resummed logarithms would certainly not be "leading". This is a consequence of the fact that QQ 
systems in the threshold regime are multi-scale problems. At the NNLO level, where all interactions can be treated as 
instantaneous, only the relative momentum of the top quarks ~ Mta s and the top quark mass Mt are relevant scales. 
At NNNLO, however, also the energy of the top quarks ~ Mta 2 a arises as a relevant scale and leads to a much more 
complicated structure of the anomalous dimensions. It is in fact not clear whether not even lower scales Mta™ (n > 2) 
become relevant if effects beyond the NNLO level are considered. 
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Figure 1: The NNLO Abelian corrections to the cross section in percent for Tt = 1.56 GeV (solid 
line) and 0.80 GeV (dashed line). See the text for more details. 



Conclusion 

Due to the large top quark mass the tt system in the threshold regime offers a unique opportunity to 
study strong interactions in heavy-quark-antiquark pairs in the threshold regime using perturbative 
QCD. In this talk I have shown that NRQCD provides an ideal framework to determine the tt cross 
section at threshold at future lepton colliders. NRQCD, an effective field theory of QCD, allows for the 
calculation of the cross section (including various distributions) to in principle arbitrary precision by 
offering a systematic formalism which parameterizes all higher order effects from first principles QCD. 
In this respect NRQCD is superior to the present day potentialmodel-like approach used for analyses 
of ti production at threshold because NRQCD does not necessarily rely on any phenomenological 
input. I therefore propose that the potentialmodel-like approach should eventually be abandoned. In 
this talk I have given a detailed recipe how NNLO relativistic corrections to the total vector current 
induced cross section can be calculated using NRQCD, and I have presented explicit results for all 
Abelian NNLO contributions. 

This work is supported in part by the U.S. Department of Energy under contract No. DOE DE-FG03- 
90ER40546. 
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